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ABSTRACT 

We discuss the deep galaxy counts from the Hubble Deep Field (HDF) imaging survey. 
At faint magnitudes, the slope of the differential number-magnitude relation is flatter 
than 0.2 in all four HDF bandpasses. In the ultraviolet, a fluctuation analysis shows 
that the flattening observed below [/300 ~ 26 mag is not due to incompleteness and is 
more pronounced than in the other bands, consistent with the idea that a redshift limit 
has been reached in the galaxy distribution. A reddening trend of s=s 0.5 magnitude is 
observed at faint fluxes in the colour-magnitude diagram, (L/300 — V606)e// versus Vqqq. 
We interpret these results as the effect of intergalactic attenuation on distant galaxies. 
At flux levels of AB sa 27 mag and in agreement with the fluctuation analysis and the 
colour-magnitude relation, about 7% of the sources in E/300, 30% in -B450 and 35% in 
V606 are Lyman-break "dropouts" , i.e. candidate star-forming galaxies at z > 2. By 
integrating the number counts to the limits of the HDF survey we find that the mean 
surface brightness of the extragalactic sky is dominated by galaxies that are relatively 
bright and are known to have (z) ~ 0.6. To AB w 29 mag, the integrated light from 
resolved galaxies in the /-band is 2.1^03 x 10 -20 erg cm -2 s _1 Hz -1 sr -1 , and its spec- 
trum is well described by a broken power-law (J„ oc A 2 from 2000 to 8000 A and I v cx A 
from 8000 to 22000 A). We discuss the predictions for the counts, colours, and lumi- 
nosity densities from standard low-go pure-luminosity-evolution models without dust 
obscuration, and find that they are unable to reproduce all the observed properties of 
faint field galaxies. 
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- quasars: absorption lines - ultraviolet: galaxies 



1 INTRODUCTION 

As the best view to date of the optical sky at faint flux lev- 
els, the Hubble Deep Field, (HDF) imaging survey has rapidly 
become a key testing ground for models of galaxy evolution. 
With its depth - reaching 5-cr limiting AB magnitudes of 
roughly 27.7, 28.6, 29.0, and 28.4 in the F300W, F450W, 
F606W, and F814W bandpasses (Williams et al. 1996) - and 
four-filter strategy in order to detect Lyman-break galax- 
ies at various redshifts, the HDF offers the opportunity to 
study the galaxy population in unprecedented detail. In 
particular, the great cosmological importance of measuring 
the number-magnitude relation (and its first moment, the 
integrated extragalactic background light) and the colour 
distribution at the faintest detection limits reached by the 
HDF stems from the prospect of using them as a probe of 



the cosmic history of star formation in galaxies (Metcalfe et 
al. 1996; Ferguson & Babul 1997). 

At faint magnitudes, the interpretation and detailed 
modeling of the observations require the self-consistent in- 
clusion of the effect of intergalactic attenuation on galaxy 
counts and colours. Absorption by intervening material 
has been known for quite some time to distort our view 
of objects at cosmological distances. It has been realized 
only recently, however, that the increasing opacity of the 
intergalactic medium (IGM) at high redshifts can be effi- 
ciently used to identify galaxies at z > 2. As shown by 
Madau (1995) (see also Yoshii & Peterson 1994), the accu- 
mulated Lyman-continuum absorption from the Lycv forest 
clouds and Lyman-limit systems along the path is so se- 
vere that galaxies beyond z ~ 3 become effectively un- 
detectable in the ultraviolet. Similarly, at z » 3.5, the 
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blanketing by discrete absorption lines in the Lyman se- 
ries is so effective that galaxies appear about 1 mag fainter 
in the blue. Ground-based observations have used colour 
techniques which are sensitive to the presence of a Lyman- 
continuum break superposed to an otherwise flat UV spec- 
trum to identify or set limits on the number of high redshift 
galaxies (Guhathakurta, Tyson, & Majewski 1990; Steidel 
& Hamilton 1992; Steidel et al. 1996a). Specific colour se- 
lection criteria which exploit the ubiquitous effect of inter- 
galactic absorption and provide a robust separation between 
high-redshift and low-redshift galaxies have been developed 
for the HDF bandpasses by Madau et al. (1996). 

The magnitude and cosmological importance of inter- 
galactic absorption on galaxy number counts can be effec- 
tively illustrated by including its effect in the standard k- 
correction term needed to translate the galaxy magnitude 
at Earth into its rest-frame value. Assuming no intrinsic 
luminosity evolution, the mean fc-correction is 



k{z) = -2.5 log 



(1 + *) 



L(u e 



L(v h s ) 
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where L(f om ) is the specific power emitted by a source at 
redshift z, i> ^ s = v ejn /(l + z), and (e~ T ) is the cosmic trans- 
mission averaged over all lines of sight. Figure 1 shows the 
average fc-correction in the £/300 and B450 bands as a func- 
tion of redshift for synthetic spectra of galaxies (based on 
Bruzual & Chariot's 1993 libraries) which well reproduce the 
colours of present-day ellipticals and late-type spirals. The 
effect at high redshifts is huge. Due to intergalactic attenua- 
tion alone, the average fc-correction in the F300W bandpass 
increases by more than 4 magnitudes between z « 2 and 
z ~ 3.5, giving origin to a "UV dropout". In the F450W 
band, the increase is only noticeable above 2 w 3, but be- 
comes very large above z ~ 4 (producing a "blue dropout"). 
Qualitatively, the inclusion of intergalactic attenuation in 
any galaxy evolution model will cause a flattening, more 
pronounced in the bluest bands, of the slope of the number- 
magnitude relation and a reddening of the ultraviolet-optical 
colours. 



In this paper we will focus on the N(m) relation in the 
£^300, -B450, Vf506, and Igi4 HDF bandpasses. A fluctuation 
analysis on the t/300 image allows us to estimate the slope 
of the ultraviolet number counts below the nominal detec- 
tion limit. These are found to flatten significantly below 26 
mag, while a reddening trend of « 0.5 mag is observed in 
the colour- magnitude diagram, (t/300 — V606)eff versus V<506. 
Both these results are consistent with the idea that, at flux 
limits fainter than 27 mag, high-z galaxies make a signifi- 
cant contribution to the counts in the redder bands, in agree- 
ment with an extension to faint magnitudes of the "dropout" 
colour technique. We compute the extragalactic background 
light (EBL) from discrete sources, and find that it can be 
well described by a broken power-law, I v oc A 2 between 2000 
and 8000 A , and I v oc A between 8000 and 22000 A. Rel- 
atively bright galaxies at the limit of existing spectroscopic 
surveys are responsible for about 60% of the sky bright- 
ness. Finally, we examine a representative (Ho — 50 km 
s- 1 Mpc- 1 , q a = 0.05) pure-luminosity-evolution (hereafter 
PLE) model without dust obscuration, and show that this is 
unable to reproduce all the observed properties of faint field 
galaxies. In particular, it overpredicts the fraction of z > 3 
objects and fails to account for the steep observed trend of 
the comoving luminosity density in the interval < z < 1 
at all wavelengths. 

Throughout this paper, all magnitudes will be given in 
the AB system (Oke 1974). 

2 THE HUBBLE DEEP FIELD 

The HDF observations allow us to statistically analyze the 
high-z galaxy population at magnitudes much fainter than 
current spectroscopic limits. In this section, we shall use 
the number counts, the Lyman-break dropout technique, 
and the colour-magnitude diagram to assess the effects of 
an absorption-induced loss of sources in the ultraviolet and 
blue bands. 

The galaxy sample we use, extracted from Version 2 
of the HDF catalog, consists of 2819 objects detected in 
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Figure 2. Differential U, B, V, and / galaxy counts as a function of AB magnitudes. The sources of the data points are indicated in each 
panel. Note the decrease of the logarithmic slope d log N/ dm at faint magnitudes. The flattening is more pronounced at the shortest 
wavelengths. The lines show the predicted counts for a PLE model (Hq = 50 Km s^ 1 Mpc - 1 , qo = 0.05), with (solid line) and without 
(dotted line) the inclusion of intergalactic absorption. The shaded region in panel (a) shows the results of our "fluctuation analysis" (see 
text). 
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the three WF chips from the F606W+F814W summed im- 
ages. The detection threshold was set to 4:a B k y , and the 
minimum area to 25 drizzled pixels (Williams et al. 1996). 
Magnitudes or upper limits in each band were computed 
from an area corresponding to the limiting isophote of the 
summed F606W+F814W image. Details of the data reduc- 
tion, source detection algorithm, and photometry are given 
in Williams et al. (1996), together with the angular separa- 
tion and colour criteria adopted for the merging of subcom- 
ponents. It is important to keep in mind that the results 
discussed below refer to a red-selected sample: in principle, 
objects with very blue colours, due for example to strong 
emission lines in the F300W or F450W bandpass, may have 
escaped detection if their red Veoe + Isi4 magnitudes do not 
satisfy the selection criteria. 

Figure 2 shows the HDF galaxy counts as a function of 
AB isophotal magnitudes in the F300W, F450W, F606W, 
and F814W bandpasses, as derived by Williams et al. (1996) 
for all galaxies with signal-to-noise ratio S/N > 3 within 
the band. At this limit the counts are likely to be more 
than 80% complete. A compilation of existing ground-based 
data is also shown. All the surveys have been corrected to 
AB magnitudes, while the second order colour corrections 
for the differences in the filter effective wavelengths have 
not been applied to the ground-based data. For the typ- 
ical colours of galaxies in the HDF, these corrections arc 
less than 0.1 mag. The HDF counts agree reasonably well 
with previous surveys, to within 20% in the magnitude range 
22 < AB < 26. At faint fluxes, however, the Williams et 
al. counts in the four individual passbands are lower by 
about 30% than those obtained for the HDF by Metcalfe 
ct al. (1996). A detailed comparison of the Williams et 
al. and Metcalfe et al. catalogs reveals several significant 
differences. First, there are systematic discrepancies in the 
total magnitudes of the galaxies. At V^oe = 25 ± 0.5, the 
median shift (Metcalfe— Williams) is —0.01 mag. This grows 
to about —0.20 mag at fainter fluxes (Viscm = 28-^29), and is 
probably due to different algorithms used for "growing" the 
photometry beyond the outer isophotes of the galaxies. A 
further difference between the two catalogs is the amount of 
splitting near bright objects, as these are split into slightly 
more pieces in the Metcalfe et al. catalog than in Williams et 
al. Finally, the Metcalfe et al. sample contains sources with 
lower S/N than the Williams et al. sample. As the error on 
the magnitudes near the detection threshold is ~ 0.5 mag, 
the slightly fainter detection limit will tend to steepen the 
counts, even a magnitude or so above the nominal depth of 
the catalog. These three effects probably account for most 
of the discrepancies, and we believe are indicative of the pos- 
sible systematic errors that are inherent in HST faint-galaxy 
photometry. 

It is important to note that, in all four HDF bands, the 
slope of the differential number-magnitude relation is flat- 
ter than 0.2 below AB = 26 mag, as previously found by 
other authors (Williams et al. 1996; Metcalfe et al. 1996), 
and that this flattening is more pronounced at the shorter 
wavelengths. We will show in Section 2.2, through a fluctu- 
ation analysis, that the turnover observed in the t/300 band 
is not due to incompleteness, but it is consistent with the 
idea that, at faint magnitudes in the red bands, a signifi- 
cant fraction of galaxies are actually at high redshifts, and 
are therefore "reddened" by absorbing material along the 



line of sight. 

2.1 Biases 

We should caution the reader that there are a number of 
subtle effects that can affect the photometry at very faint 
magnitudes, especially in the f/300 band. First, the charge 
transfer efficiency (CTE) of the WFPC-2 CCD is not per- 
fect, and can lead to magnitudes that are systematically 
too faint for sources observed at low count rates with low 
background. For the typical background levels of 1 DN for 
the F300W exposures, the average CTE correction is ~ 0.03 
magnitudes for galaxies with U300 ~ 25. This correction 
will depend on the location of the source on the chip, vary- 
ing from zero at low row numbers to 0.5 mag near the top 
of the chip. The CTE effect on sources several magnitudes 
fainter is not well calibrated, but is unlikely to be more than 
0.1 mag. A second subtlety is that, at background levels of 
~ 1 DN per exposure, the determination of the local back- 
ground for a galaxy is affected by the discrete steps of the 
CCD analog-to-digital converter. While the algorithm used 
by FOCAS (a clipped mean) should be a good estimator of 
the background, we have not tested this ourselves on simu- 
lated data. Finally, for galaxy fluxes within several percent 
of the background, it is possible that systematic biases creep 
in due to the different treatment of the pixel intensities in the 
source aperture and background region. In particular, the 
mean intensity in the local-background region is determined 
after clipping out pixels more than 10 sigrna away from an 
initial estimate of the mean, while the flux in the source 
aperture is an undipped sum of the counts. This could bias 
the fluxes of sources that are very faint in F300W towards 
slightly higher values, resulting in measured colours some- 
what bluer than the true colours. Finally, at magnitude 
levels fainter than Vooa = 28, it is likely that even objects 
which are brighter than the nominal depth of the catalog 
may be missed due to low surface brightness. These galax- 
ies could be preferentially either blue or red. For example, 
high redshift galaxies that are not forming stars rapidly will 
be red and have quite low surface brightness. The (1 + z) 4 
dimming effect is equal to 1.2 magnitudes between z = 2 
and z — 3. On the other hand, locally, LSB galaxies tend to 
be quite blue. A relatively local population of LSB galax- 
ies could in principle exist and the fainter examples could 
disappear below the HDF sky background. 

With these caveats in mind, we will limit our analysis 
to galaxies with Veoe < 29 mag, with the expectation that 
systematic biases will be significantly less than ~ 0.5 mag in 
the photometry, and less than ~ 50% in the galaxy counts. 

2.2 Fluctuation Analysis 

As the F300W bandpass is sensitive to intergalactic atten- 
uation at z > 2, deep counts in this band are essential for 
estimating, through the comparison with counts at longer 
wavelengths, the number of high- 2 galaxies in a red-selected 
sample, and for deriving constraints on models of galaxy 
evolution. 

The published UV N(m) (Williams et al. 1996) al- 
ready suggests a sharp flattening from a logarithmic slope 
of dlog N/dm « 0.40 at f/300 < 26 mag to a slope consistent 
with zero for 26 < [/300 < 28. In principle, however, this 
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Figure 3. S/N distribution in the t/300 band. The data (shaded 
region) are compared to a simulated distribution with U c = 25.75 
mag, and ■yf = 0.135 (dashed line) or 0.350 (dotted line). The 
gaussian curve represents the expected distribution for unde- 
tected objects (see text). The upper right panel shows the 95% 
confidence region for the parameters U c and ■yf from a KS test 
(see text). 



sharp flattening could also be due to significant incomplete- 
ness in the faintest C/300 bins. Note that, since the limiting 
ultraviolet flux is about one magnitude brighter than in the 
other bands, the fraction of galaxies in the catalog that are 
undetected (i.e. S/N < 3) in the F300W bandpass is very 
large (~ 80%), and it is difficult to draw statistically reliable 
conclusions from the raw data. 

In the attempt to push the I/300 counts to even fainter 
magnitudes, we have therefore applied a "fluctuation anal- 
ysis" to the Version 2 catalog data. The simplest possi- 
ble model for the differential counts has been assumed, i.e. 
a double power-law characterized by four parameters: the 
overall normalization (A), two slopes (75 and 7/ at bright 
and faint magnitudes, respectively), and the magnitude (U c ) 
of the turnover. While two of them (A and 75) have been 
kept fixed by requiring a good fit to the bright HDF N(m), 
we have determined a best estimate for the other two pa- 
rameters (7/ and Uc) in the following way: 

a) values for the ultraviolet magnitudes of the sources 
have been randomly extracted from the input N(m); 

b) the corresponding "measured" counts for each galaxy 
have been derived assuming Poisson statistics; 

c) an isophotal area has been assigned to each object, 
on the basis of the correlation between area and magnitude 
observed in the catalog, and the corresponding number of 
background counts has been derived; 

d) finally, a signal-to-noise ratio [S/N = Li /a (Li)] has 
been computed according to the formula (Note that in equa- 
tion 6 in Williams et al. 1996 there was a typo in the second 



and third terms): 

[Vo(Li)] 2 =TN obj + 1.9 2 T 2 cr 2 3ky A obj + 

2 2 2 2/ V / 

1.9 r CT s ky A ffj / A s ky , 

where Li is the sky-subtracted number of counts within the 
detection isophote, F is the inverse gain, N tj is the total 
number of counts in the object aperture, A ob j and A sky are 
the areas within the source and sky apertures, and a s k y is 
the measured standard deviation of the background within 
the sky aperture (cf. Williams et al. 1996). 

By repeating these steps a large number of times 
(10,000), we have constructed a simulated distribution of 
signal-to-noise ratios for each input N(m), and these distri- 
butions have then been compared to the observed one for all 
the 2819 objects of the merged catalog. Figure 3 shows such 
a comparison for U c = 25.75. As expected, the observed dis- 
tribution can be described as a gaussian at low S/N values, 
which represents the expected distribution for undetected 
objects, plus a tail at high S/N values. Since these are com- 
puted not in random positions, but at the positions where 
galaxies are detected in the red bands, the gaussian-like 
distribution is not centered around zero, but is displaced to- 
ward a higher value (S/N ~ 0.5). The two curves plotted 
in Figure 3 show the results of two simulations, correspond- 
ing to the slopes of 7/ = 0.380 and 7/ = 0.135. The case 
7/ ~ "/b ~ 0.4 is clearly ruled out to a high degree of signifi- 
cance. From a KS statistical test the formal 95% confidence 
region for (£4,7/) can actually be derived. The two param- 
eters are obviously highly correlated, with steeper slopes 
allowed for brighter values of U c and vice- versa (see inset in 
the upper right corner of Figure 3). Note, however, that the 
region in the upper left part of the 95% confidence ellipse 
produces number counts which are lower than the observed 
ones for 25.75 < C/300 < 27.25, while its lower right part pro- 
duces counts which are higher than the observed ones in the 
same magnitude range. By requiring the models to be con- 
sistent with the data to better than 5% in the four bins be- 
tween 25.75 < t/300 < 27.25, the two hatched regions of the 
parameter space can be further excluded, and the allowed 
ranges become 25.6 < U c < 26.2 and 0.08 < 7/ < 0.16. 
The best-fit values are U c = 25.75, and 7/ = 0.135. The 
latter is steeper than the value of 0.05 quoted in Williams 
et al. (1996), but still flatter than the corresponding values 
for the -B450, Veoe, L314 bands, as expected in the case of a 
significant absorption-induced loss of faint, distant sources. 

All the allowed models predict a number of galaxies be- 
tween 20 to 30% higher than observed in the faintest bin 
quoted by Williams et al. (1996), 27.5 < t/300 < 28.0. The 
Version 2 catalog contains 112 galaxies in this magnitude 
bin, with only 79 satisfying the S/N > 3 constraint used 
by Williams et al. in deriving their counts. This suggests 
some incompleteness in the Williams et al. counts at this 
magnitude level. Since the catalog has been derived using 
isophotal magnitudes and the isophotal areas are defined in 
the red bandpasses, for galaxies with a given ultraviolet flux 
(i.e., with a given signal) there is an inverse correlation be- 
tween S/N and isophotal area: smaller signal-to-noise ratios 
correspond to larger areas and therefore, in first approxima- 
tion, to galaxies that are brighter in the red bands. For 
this reason, the sources eliminated because of the S/N con- 
straint are not a random sub-sample of the galaxies in this 
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magnitude range, but correspond to those with the reddest 
C/300 — V(jo6 and/or C/300 — Isn colours. This has to be taken 
into account if one analyzes the colour distribution of the 
galaxies in the faintest C/300 magnitude bin of the catalog. 

As mentioned in Section 2, galaxies with very blue 
colours, not selected in the catalog because fainter than the 
limits in V6O6 + ^814, could in principle contribute some- 
what to the C/300 counts. Because our analysis uses as input 
the red selected catalog, this contribution would not be ac- 
counted for in our extrapolated estimate of the t/300 counts. 
Given the approximate limit in the Vk>6 band (Vn m ~ 29.0), 
these "potential" blue objects could appear at t/300 ~ 29 if 
they had C/300 — Veoe < and at C/300 ~ 29.5 if they had 
C/300 — Veo6 < 0.5. In the "bright" part of the catalog (i.e. 
V606 < 26), where no such object is expected to have been 
lost, there is only 1% of the galaxies with C/300 — Veoe < 
and only 12% of the galaxies with C/300 — Vfjo6 < 0.5. On 
this basis we argue that our conclusion on the reality of the 
flattening in the C/300 counts is robust and our estimate of 
the parameters of the slope down to C/300 ~ 29.5 can be 
considered to be reliable. 

2.3 Lyman-Break Dropouts and the High-z 
Luminosity Functions 

Photometric criteria for robustly selecting Lyman-break 
galaxies in the HDF have been developed by Madau et al. 
(1996). They have been tuned up to provide what appear 
to be largely uncontaminated samples of star-forming galax- 
ies at high redshifts, and have been refined (Madau 1997) 
after the many redshift measurements with the Keck tele- 
scope (Steidel et al. 1996b; Cohen et al. 1996; Lowenthal 
et al. 1997). The HDF ultraviolet passband - which is 
bluer than the standard ground-based U filter - permits 
the identification of star-forming galaxies in the interval 
2 < z < 3.5. Galaxies in this redshift range predominantly 
occupy the top left portion of the C/300 — -B450 vs. B450 — isi4 
colour-colour diagram because of the attenuation by the in- 
tergalactic medium and intrinsic absorption. Galaxies at 
lower redshift can have similar C/300 — -B450 colours, but 
are typically either old or dusty, and are therefore red in 
B450 — ^814 as well. In analogous way, it is possible to single 
out star-forming galaxies at 3.5 < z < 4.5 by looking at the 
-B450 — V<306 vs. V606 — Csi4 colour-colour diagram. 

We have estimated the number of C/300 and B450 
dropouts as a function of V6O6 magnitudes by applying, with 
some modifications, the procedure described in Madau et al. 
(1996). The main differences with respect to Madau's pro- 
cedure are: i) We considered as detected only objects with 
S/N > 2.5, rather than S/N > 1 as in Madau's analysis. 
For the objects with a lower S/N we computed a lower limit 
to the magnitude corresponding to S/N — 2.5. As seen in 
Figure 3, it is only for S/N > 2.5 that the observed dis- 
tribution starts to depart significantly from the expected 
distribution of the noise, represented in the figure by the 
gaussian curve. Therefore, the objects with lower S/N can 
not be considered detections, ii) Since with this choice we 
have a higher number of lower limits in the colours, the num- 
ber of dropouts in each V6O6 bin has been estimated using 
the Kaplan-Meier estimator in the package ASURV (La Val- 
ley, Isobe & Feigelson 1992), which implements the methods 
for a maximum-likelihood reconstruction of the true distri- 
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Figure 4. Fraction of high-z galaxies as a function of V606 
magnitudes as obtained from a dropouts analysis, using the 
ASURV package. Filled circles: C/300-dropouts. Empty circles: 
C?450-dropouts. Solid lines: PLE model predictions. 

bution function when limits are present (Feigelson & Nelson 
1985). Table 1 shows the numbers of C/300 dropouts for 
Veoe < 27.5 and of B450 dropouts for Veoe < 28 as a func- 
tion of V606 magnitudes. At fainter V6O6 magnitudes the 
number of limits in C/300 becomes greater than 60% of the 
total number of data points. For these high percentages of 
limits the performances of survival analysis statistics dete- 
riorate and the results are not reliable anymore (Feigelson 
1992). The corresponding fractions of dropouts are shown 
in Figure 4 as a function of Veoe magnitudes in 0.5 magni- 
tude bins. The fraction of U + B dropouts varies from ~ 5% 
to 40% in the range 23.5 < Veo6 < 27.5. The cumulative 
fraction of F300W-dropouts to Veoe < 27.5 is 28 ±2%, while 
it is 2.5 ± 0.6% to Veoe < 28 for the F450W-dropouts. 

With a similar procedure we have computed the fraction 
of Lyman-break dropouts as a function of C/300 and B450 
magnitudes, finding that this decreases significantly in the 
C/300 band. For example, at AB m 27, it is ~ 7% in C/300, 
while it is ~ 30% in B450 and ~ 35% in Veoe- 

The Kaplan-Meier estimator used in this analysis can 
be properly applied only in presence of a random censorship, 
i.e. in cases in which the distribution of the observational 
sensitivities for the non-detected objects is the same as that 
for the detections (Schmitt 1985, Magri et al. 1988). In our 
case this condition is satisfied, because for each magnitude 
bin in V6O6 the distribution of the thresholds in C/300 is es- 
sentially a function of A ob j (see Eq. 2) , which is not related 
to the variable to be measured (i.e. the C/300 magnitude). 
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However, the non-parametric Kaplan-Meier reconstruc- 
tion of the true distribution implicitely assumes also that the 
limits and the detections are drawn from the same parent 
distribution. Since we can not be sure that this assump- 
tion is true, we have analyzed the same data also with a 
maximum likelihood parametric method. Consistently with 
the analysis described in more detail in Section 2.4, we have 
assumed the following functional form for the distribution 
of the C/300 — -B450 colours: a fraction / of the population 
of galaxies is highly absorbed in C/300 and is described for- 
mally by a (5-function at a high value of C/300 — -B450 ; a frac- 
tion (1 — /) of the population is described by a gaussian 
distribution in C/300 ~ #450- The likelihood for this func- 
tional form can be easily derived (see Avni & Tananbaum 
1986). Maximizing the likelihood function one can then ob- 
tain the best estimates for the three parameters (/ plus the 
two describing the gaussian) and from these derive the best 
estimate for the number of dropouts, i.e. those satisfying 
the condition of having a C/300 — S450 color greater than the 
adopted threshold. The numbers derived with this proce- 
dure are all consistent, within one sigma, with those given 
in Table 1 in each magnitude bin. Since the assumptions 
in the two methods are significantly different, this suggests 
that our estimate of the number of dropouts is statistically 
robust. 

As a test of the adopted colour selection window, we 
have compared our estimates of the fraction of bright UV 
dropouts in the HDF with the results from the B-selected 
spectroscopic sample described by Cowie et al. (1997). The 
latter covers, in the magnitude range 24 < B < 24.5, about 
13 arcmin 2 , and contains 42 galaxies, of which 4 have 2 > 
2, and 2 are still unidentified (and most likely lies at z > 
1.6, Cowie et al. 1997). The fraction of spectroscopically 
identified galaxies at z > 2 in the ground-based sample is 
therefore between 10 and 14%. This is in good agreement 
with our estimated fraction, ~ 11%, of ultraviolet dropouts 
in the same magnitude range. 

Using the numbers of dropouts as a function of mag- 
nitude, we have computed the luminosity functions (LF) of 
high-2 galaxies in two different redshift ranges. The ob- 
served Vf506 and 7gi4 magnitudes have been converted to 
intrinsic M1600 AB magnitudes, and the number densities 
of C/300 and B450 dropouts have been obtained assuming 
that they sample the redshift intervals 2.0 < 2 < 3.5 and 
3.5 < z < 4.5, respectively. The results are shown in Figure 
5 {H = 50 km s _1 Mpc -1 , q = 0.05). The shape of the 
LF for the redshift interval 2.0 < 2 < 3.5 is reasonably well 
defined and can be described by a Schechter function with 
parameters a = -1.3, M* 600 = -21.7 and $* = 0.8 x 10~ 3 
Mpc -3 . For the higher redshift interval the fitting line 
shown in the figure has the same a and M* 600 parameters, 
but a lower normalization ($* = 0.2 x 10~ 3 Mpc -3 ). As 
a comparison, the thick curve in Figure 5 shows the local 
M1600 LF, derived from the Treyer et al. (1997) data. The 
high-z luminosity functions clearly extend to much brighter 
magnitudes, while no indication of a steep faint end slope is 
clearly seen in the HDF data. Even the HDF data, however, 
are not deep enough to exclude the possibility of a steepen- 
ing at faint fluxes. The comoving luminosity density result- 
ing from the integration to M1600 = —18 of the two HDF LF 
are p L = 2.0(±0.5) x 10 19 WHz -1 Mpc -3 for 2.0 < z < 3.5, 
and p L = 5.2(±2.6) x 10 18 WHz -1 Mpc -3 for 3.5 < z < 4.5, 
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Figure 5. Luminosity function at 1600 A (Ho = 50 km s _1 
Mpc -1 , go = 0.05). Filled circles: C/300-dropouts. Empty circles: 
i?450-dropouts. The solid lines represent Schechter fits to the 
data. The local luminosity function at 2000 A from Treyer et al. 
(1997) is also shown (thick curve). 

in good agreement with the estimates of Madau (1997). The 
errors on pl have been obtained by combining the statisti- 
cal errors on the number of dropouts (see Table 1) with an 
estimated uncertainty of the volumes effectively sampled by 
the adopted colour selection techniques. 

These luminosity functions (and the corresponding lu- 
minosity densities) should probably be considered as lower 
limits to the true ones. In fact, on the basis of an analysis of 
the spectra and colours of a representative, spectroscopically 
confirmed sample of high redshift star-forming galaxies, Pet- 
tini et al. (1997) have shown that the likely dust correction 
to the integrated ultraviolet luminosity for these galaxies is 
about a factor of 3. (Note that the Treyer et al. LF has not 
been corrected for extinction). Moreover, the dropout tech- 
nique used in selecting our high redshift candidates works 
efficiently for galaxies for which the attenuation due to in- 
tergalactic matter is equal to or higher than the assumed 
average attenuation. Galaxies with a lower than average at- 
tenuation may well have (C/300 — -B450) and (-B450 — V&o&) 
colours bluer than those used in selecting the high-2 sam- 
ple. Madau et al. (1996) have estimated that only about 
20% of galaxies at 2 > 2.7 can be missed because of this 
effect. However, only extensive spectroscopic observations 
of galaxies lying below the colour threshold can provide a 
quantitative estimate of the size of this incompleteness. 

2.4 Colours 

As already mentioned, the number of objects detected in 
C/300 is a relatively small fraction of the total number of 
galaxies in the catalog. This obviously prevents a detailed 
analysis of the full C/300 — V606 distribution at faint Veoe 
magnitudes, where for most of the objects only an upper 
limit to the ultraviolet flux is available. 
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Figure 6. Colour-magnitude diagram, ((7 — V) e ff vs. V606 (filled 
circles). The six brighter bins contain 80 galaxies each, while the 
fainter bins contains 240 objects. The dashed curve represents a 
quadratic fit to the data to Veoe < 29. Solid line: prediction from 
a PLE model including the attenuation by the IGM. Dotted line: 
unattcnuated PLE model. 

It is possible, however, to use the measured fluxes after 
background subtraction for all the galaxies to get a statis- 
tical estimate on the average properties of the C/300 — Vfjoe 
vs. Vfjoe relation at much fainter V<506 magnitudes than al- 
lowed by the C/300 detections. We have therefore computed 
in each bin of Vf306 magnitude the total F300W and F606W 
fluxes as the sum of all the counts measured in individual 
galaxies, and have derived an "effective" colour for the en- 
tire population, (C/300 — V6oe)eff. This is not equal to the 
mean value, (C/300 — Hoe)- For example, if the colour dis- 
tribution were a gaussian with mean (C/300 — Veoe) and dis- 
persion a, it is easy to derive the relation (C/300 — Vfjoe) = 
(C/300 - Veoe )eff + 0.4605a 2 . 

Figure 6 shows the resulting colour-magnitude diagram. 
After an initial blueing trend, which is expected in PLE 
models (Pozzetti, Bruzual & Zamorani 1996), at fainter 
magnitudes galaxies become redder in (C/300 — Vfjo6)eff • This 
may be due to an increase in the fraction of high redshift 
galaxies or to a change in the colour properties of low red- 
shift objects. As a test for distinguishing between these two 
alternative possibilities we have estimated the expected ul- 
traviolet counts from the observed optical counts and the 
derived (C/300 — V606)eff colours. The following input ingre- 
dients have been used: 

i) the observed log N(m) — Vqo6 relation, which we fitted 
as a double power-law with 71, = 0.32, 7/ = 0.17, and V c = 
25.5; 

ii) an assumed gaussian distribution for C/300 — V606 with 
a mean such that the resulting (C/300 — V606)eff is equal to 
the fitted value in each V606 bin; 

iii) two different assumptions for the fraction, /(Vtsoe), 
of high redshift galaxies in 0.5 magnitude bins: /i(Vf506) = 
and / 2 (VWs) = 0.097 x (Veoe - 22.68) for 23 < Veoe < 27, as 
obtained by the dropout technique (see § 2.2). At fainter flux 
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Figure 7. C/300 (open triangles) and V6O6 (filled triangles) num- 
ber-magnitude relation from the HDF data. Simulated ultraviolet 
counts arc also plotted, with (open circles) and without (filled cir- 
cles) the inclusion of high-z galaxies (see text). The solid lines 
depict the best-fits to the Veoe data and the t/300 N(m) relation, 
with the former brightened by 1 magnitude for clarity. 

levels, where the dropout analysis fails, we have assumed 
that /2(V606) remains constant. 

Figure 7 shows the results for the two cases. In the 
very simple approach followed here we have assigned the 
same colour, C/300 — Veoe = 4 mag, to all high redshift ob- 
jects, which therefore contribute only a negligible fraction 
of the total measured C/300 flux in each V606 bin. The 
(C/300 — V606) e ff colour for the low redshift galaxies is then 
recomputed by scaling the total measured V606 flux in each 
bin by 1 — /2(V606). In the same figure we show the data 
points with error bars, the fit to the bright counts, and the fit 
at fainter magnitudes derived from the fluctuation analysis. 
The predicted ultraviolet N(m) is significantly higher than 
what is allowed by the fluctuation analysis for / = /1 , while 
it is in excellent agreement with the fluctuation analysis for 
/ = MVeoa). 

We conclude that the results of our analyses, i.e. the 
fluctuation analysis on the C/300 data, the colour selection 
of Lyman-break galaxies, the derivation of the C/300 — V506 
effective colour as a function of Veoe magnitude (and its 
application to the V6O6 counts to predict the expected ultra- 
violet counts), are all consistent with each other and with 
the idea that, at faint UV flux limits, there is a significant 
absorption-induced loss of sources in the HDF. 

3 EXTRAGALACTIC BACKGROUND LIGHT 

The extragalactic background light (EBL) is an indicator 
of the total optical luminosity of the universe. It can pro- 
vide unique information on the origin of structures at early 
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Table 2. Integrated Galaxy Light 
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Figure 8. Extragalactic background light per magnitude bin as a 
function of U (filled circles), B (open circles), V (filled triangles), 
I (open squares) and K (filled squares) magnitudes. For clarity, 
the B, V, I and K values have been multiplied by a factor of 2, 
8, 30, and 100, respectively. 

epochs, as the cumulative emission from pregalactic, pro- 
togalactic, and evolving galactic systems is expected to be 
recorded in this background (Tyson 1995). The contribu- 
tion to the EBL from discrete objects can be calculated 
directly by integrating the flux times the differential num- 
ber counts down to the detection threshold. We have used 
the compilation of ground-based and HDF data shown in 
Figure 2, and the compilation of 7f-band data shown in 
Figure 3 of Pozzetti et al. (1996), to compute the EBL at 
3500 < A < 22000A. The results are listed in Table 2, along 
with the magnitude range of integration and the estimated 
la error bars, which arise mostly from field-to-field varia- 
tions in the numbers of relatively bright galaxies. 

Because of the flattening at faint magnitudes of the 
N(m) differential counts, most of the contribution to the 
EBL comes from relatively bright galaxies. This is clearly 
seen in Figure 8, where the function i„ = io _0 - 4 ( m + 48 - 6 ) x 
N(m) is plotted against apparent magnitude in all bands. 
While the differential contribution to the EBL increases at 
bright fluxes, where the slope of the ground-based counts is 
relatively steep, it appears that the HDF survey, reaching 
the limits beyond which the counts flatten below 0.40, has 
achieved the sensitivity to capture the bulk of the extra- 
galactic light from discrete sources. An extrapolation of the 
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Figure 9. Spectrum of the extragalactic background light as de- 
rived from a compilation of ground-based and space-based galaxy 
counts. Also plotted are the 90% all-sky-photometry upper limits 
from Bowycr (1991), Toller (1983), Dube et al. (1977, 1979), and 
Boughn et al. (1986). For clarity, the upper limit at A = 2.2fim 
has been plotted at 2Afim. Dotted line: broken power-law fit to 
the EBL spectrum. 

observed N(m) to brighter and/or fainter magnitudes would 
typically increase the integrated light by less than 20%. The 
spectrum of the EBL (including a UV point at 2000 A from 
Armand, Milliard, & Deharveng 1994) is shown in Figure 
9, and can be well described by a broken power-law: I v — 
1.9 x 10~ 2() (A/8000 A) 2 erg cm^s" 1 Hz^sr" 1 from 2000 to 
8000 A, and 7^ = 1.9x 10~ 20 (A/8000 A) erg cm" 2 s" 1 Hz" 1 sr" 
from 8000 to 22000 A. Several diffuse EBL upper limits - 
set by all-sky photometry - are also plotted. While in the 
optical/UV these are from 3 to 5 times higher than the con- 
tribution from known galaxies, the upper limit to the total 
diffuse if-band EBL lies only slightly above the summed 
flux from discrete objects. 



4 COMPARISON BETWEEN THE DATA AND 
A STANDARD PLE MODEL 

Qualitatively, the effects of the inclusion of intergalactic at- 
tenuation in any galaxy evolution model are easily described. 
By severely obscuring sources above z w 2 in the ultraviolet 
band, and above 2 ~ 4 in the blue, the ubiquitous presence 
of H 1 along the line of sight to distant sources gives effec- 
tively origin to galaxy samples that are volume-limited in 
these bands. This truncates the galaxy redshift distribution 
and causes a flattening, more pronounced in the bluest fil- 
ters, of the slope of the number-magnitude relation. The 
differential effect in the various bandpasses could then, at 
least in principle, allow us to discriminate between different 
models. Indeed, a more pronounced flattening of the ultra- 
violet N(m) relative to the counts at longer wavelengths is 
likely to flag the presence of a significant population of faint 
high-z objects. Moreover, objects close to the redshift limit 
set by the opacity of the intergalactic medium should appear 
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Figure 10. Integrated luminosity per unit comoving volume as 
a function of redshift (H = 50 km s _1 Mpc -1 , qo = 0.05). The 
sources of the data points are indicated in each panel. The curves 
show the luminosity density predicted by a PLE model at 4400 
A (solid line), 2800 A (dotted line), and 1600 A (dashed line). 



quite redder. 

In the following, we shall briefly examine the predictions 
of a standard pure-luminosity-evolution model, modified by 
the inclusion of the effects of intergalactic attenuation, and 
compare them with the results obtained from the HDF data. 
Such a model, which assumes a single redshift of forma- 
tion (zf) for all galaxies and a open universe, was shown by 
Pozzetti et al. (1996) to be able to reproduce, with differ- 
ent degrees of accuracy, most of the properties of faint field 
galaxies as derived from ground-based surveys. The differ- 
ential galaxy counts as a functions of AB magnitude, as ex- 
pected within a PLE model with Zf = 6.3, are shown in Fig- 
ure 2. Intervening absorption begins to affect the ultraviolet 
counts at C/300 ^ 25 mag, where the predicted slope flattens 
to 7 < 0.3. By C/300 ~ 29, the PLE model including inter- 
galactic attenuation has 7 ~ and predicts about 2 times 
less galaxies than expected in a transparent universe. The 
median redshift for the global galaxy distribution is z med ~ 2 
at (C/300, B450, Veoe, hu) ^ (29, 28, 27, 26). At interme- 
diate magnitudes (20 < AB < 25), luminosity evolution 
causes a blueing trend in the mean colours (Pozzetti et al. 
1996. Only at fainter flux levels (> 25 mag) the effect of the 
IGM becomes significant, producing a clear trend towards 
redder ultraviolet-optical colours with increasing magnitude 
(Figure 6). 

Overall, the adopted PLE model appears to provide a 
rather acceptable fit to both the ground-based and HDF 
galaxy counts (Figure 2), as well as to the EBL spectrum 
(Figure 9). It produces, in good agreement with the ob- 
servations, a flattening at faint magnitudes which is more 
significant at shorter wavelengths. The predicted turn-over 
in the C/300 bandpass is, however, sharper than allowed by 



the fluctuation analysis, due to the presence of too many 
high-z galaxies in the model relative to data. While, as 
shown in Figure 6, the inferred reddening in (C/300 — Veo6)eS 
at magnitudes fainter than Veoe = 26 is consistent with what 
expected in a PLE model (note how the same model with- 
out attenuation predicts bluer colours than observed in this 
magnitude range), at brighter flux levels (23 < V@oe < 26) 
the synthetic colours appear to be too blue. The reason for 
this difference between model and data is not clear. 

In Figure 4, the predicted fractions of colour-selected 
C/300 and -B450 dropouts versus Veo6 magnitude are compared 
with the HDF observations. While the relative numbers of 
ultraviolet dropouts are consistent with the PLE prediction, 
there is a clear deficit of observed B450 dropouts relative to 
the model. As already pointed out by Ferguson & Babul 
(1997), it appears that a major problem of a PLE scenario 
is the large number of predicted galaxies at high redshifts. 

This is confirmed by Figure 10, which shows a compari- 
son between the observed luminosity densities as a function 
of redshift at different wavelengths and the predictions of 
the PLE model. At 4400 A the model generates a lumi- 
nosity density which increases only slightly with redshift for 
z < 1. The local comoving emissivities derived from the data 
by different authors show a large spread, of order a factor 
two, although the most recent observations (Ellis et al. 1996; 
Zucca et al. 1997) suggest a higher normalization than that 
obtained by Loveday et al. (1992). If this is correct, the 
data are consistent with no or weak evolution up to z ~ 0.3. 
A rapid evolution, not reproduced by the model, is instead 
required by the data in the range 0.3 < 2 < 1.0, both in the 
blue and at ultraviolet wavelengths. At higher redshifts, the 
model is in good agreement with the luminosity density we 
have derived at 1600 A in the range 2.0 < 2 < 3.5, while 
it is higher than the observations by a factor ~ 10 in the 
interval 3.5 < 2 < 4.5. 



5 CONCLUSIONS 

By applying a simple fluctuation analysis to the C/300 data, 
we have obtained an estimate of the slope of the ultraviolet 
counts, down to more than a magnitude fainter than the 
last data point published by Williams et al. (1996). The 
derived UV number-magnitude relation agrees with what 
has been obtained by convolving the observed Veoe counts 
with the derived (C/300 — Vaoe)eff vs. Veo6 relation, and 
assuming a fraction of high-2 galaxies consistent with the 
results of the dropout colour technique. All these results 
convincingly show that the effects of an IGM absorption- 
induced loss of sources are significant at faint ultraviolet 
magnitudes (AB > 27). At flux levels of AB « 27 mag, 
about 7% of the sources in C/300, 30% in B450, and 35% 
in V606 are Lyman-break "dropouts", i.e., candidate star- 
forming galaxies at 2 > 2. 

The rest-frame UV luminosity densities we derive at 
high redshifts using a maximum likelihood reconstruction 
of the true colour distribution are in good agreement with 
previous estimates (Madau et al. 1996; Madau 1997). We 
note that the existing data do not necessarily require a sharp 
peak of the ultraviolet luminosity density at 2 ~ 1.5, but 
may also be consistent with an approximately constant value 
over a redshift range as large as 0.8 < 2 < 3.0. 
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From our dropout sample we have also derived the 
M1600 luminosity functions in the redshift intervals 2.0 < 
2 < 3.5 and 3.5 < 2 < 4.5. The LF is reasonably 
well defined in the lower redshift interval, and can be de- 
scribed by a Schechter function with parameters a = —1.3, 
M* 600 = -21.7 and $* = 0.8 x 10" 3 Mpc" 3 . The same a 
and M* 600 values, but with a normalization four times lower, 
provide a good fit to the blue dropouts at higher redshifts. 
Both functions extend to brighter magnitudes than the local 
2000 A LF recently derived by Treyer et al. (1997). 

A detailed comparison between a simple PLE galaxy 
evolution model and the available ground-based and HDF 
data shows that the model fails to reproduce all the observed 
properties of faint field galaxies. While able to fit reasonably 
well the observed number counts, colours, and the fraction 
of ultraviolet dropouts at {2) ~ 2.75, this PLE model over- 
predicts by about one order of magnitude the number of 
blue dropouts at 2 > 3.5 and the corresponding rest-frame 
UV luminosity density. Moreover, it fails to account for the 
steep observed trend of the galaxy emissivity observed in 
the redshift range 0.3 < 2 < 1 at all wavelengths. The mod- 
eling of the emission history of the galaxy population at late 
epochs is outside the scope of this paper; we note here that 
Madau, Pozzetti & Dickinson (1997) have recently shown 
that the rapid variation of the luminosity density at small 
redshifts can be well fit by a simple stellar evolution model, 
defined by time-dependent star formation rate per unit co- 
moving volume and a universal initial mass function. At 
high redshifts, since in our PLE model a significant frac- 
tion of the emitted light is due to elliptical galaxies just 
after their first massive starburst phase, a possible solution 
to the observed discrepancy could be obtained by assuming 
that a significant amount of dust obscures the early evolu- 
tion of these massive systems. A similar conclusion has been 
reached by Franceschini et al (1997) in their analysis of a It- 
band selected sample of early-type galaxies in the HDF. Ap- 
proximate consistency with the (2) ~ 4 observations would 
require a dust extinction of E(B - V) ~ 0.2-0.23 for SMC- 
type dust and the empirical extinction-law of Calzetti, Kin- 
ney & Storchi-Bergman (1994), respectively, corresponding 
to an ultraviolet extinction of about a factor of 10. The 
current estimates of the ultraviolet extinction for the star- 
forming galaxies at 2 ~ 3 are still very uncertain and range 
from about 2 (Pettini et al. 1997) to more than 10 (Meurer 
et al. 1997). If such a large amount of dust, consistent with 
the observations of local starburst galaxies (Calzetti 1997), 
was present in star-forming objects at high-2, it would ab- 
sorb a large fraction of the optical-UV light and re-emit 
it in the far infrared band (see, for example, the models 
described in detail in Mazzei, De Zotti & Xu 1994). In 
this case, the amount of star formation derived from our 
and similar dropout analyses should only be considered as 
a lower limit. Preliminary observations of the HDF with 
ISO (Rowan-Robinson et al. 1997) have revealed about a 
dozen far-IR (6.7 and 15/xm) sources, identified with opti- 
cal galaxies. By fitting the spectra of these galaxies from 
the ultraviolet to the far-IR, Rowan- Robinson et al. (1997) 
have derived star formation rates significantly higher than 
those derived from the UV fluxes detected with HST. On 
the other hand, Madau et al. (1997) show that a model with 
a large amount of hidden star formation at early epochs will 
overproduce the cosmic metallicity at high-2 as sampled by 



the damped Lyman-cv systems (Pettini et al. 1997) . The re- 
sults discussed above demonstrate the great potential of IR 
observations for complementing the information provided by 
the rest-frame ultraviolet and for helping our understanding 
of the nature of high-2 galaxies. 
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